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We study the electromagnetic scattering by multilayered 
biperiodic aggregates of dielectric layers and gratings of con- 
ducting plates. We show that the characteristic lengths of 
such structures provide a good control of absorption bands. 
The influence of the physical parameters of the problem (sizes, 
impedances) is discussed. 



I. INTRODUCTION 

Electromagnetic absorbers and frequency selective sur- 
faces (FFS for short) have recently received an increasing 
interest. There is a growing need for electromagnetic ab- 
sorbers, and in particular for lighter, thinner and more 
highly absorbing materials. Frequency selective surfaces 
are generally made of planar screens with periodic or 
biperiodic metallizations. One generally considers two 
types of FFS: capacitive FFS are transparent at low fre- 
quencies; inductive FFS are reflecting ones. Their behav- 
ior at the resonance frequency is complementary. Capac- 
itive FFS consist of arrays of metal patches embedded 
in a dielectric structure, which may be a stratified one. 
The dielectric structure provides the mechanical support 
of the FFS. Inductive FFS consist of perforated screens. 

Such frequency selective surfaces have been considered 
by several authors p])p|Jl3[ who have proposed various ap- 
proaches for the numerical resolution of the correspond- 
ing scattering problem. Efficient numerical methods are 
now available for the analysis and design of FFS, as we 
shall show. 

The purpose of this paper is to show that the ab- 
sorption bands of such structures may be controlled by 
combining the performances of capacitive FFS and elec- 
tromagnetic absorbers. To vary the frequency response 
of a FFS, the standard method consists in varying the 
geometry of the array elements. We give efficient com- 
putational methods for analyzing this kind of structure. 
The representation of the transmitted and reflected fields 
is obtained by applying resistive boundary conditions to 
include a general surface impedance in the problem for- 
mulation. 



We provide examples of such periodic or biperiodic 
structures, whose absorption bands can easily be con- 
trolled by varying some of the characteristic lengths of 
the system. More precisely, we consider multilayers made 
of dielectric stacks and surface gratings with various 
shapes and sizes. We show that such structures yield 
absorption bands, and that the location and bandwidth 
of such bands may be controlled by varying the charac- 
teristic sizes of the structure. 

This paper is organized as follows. After this introduc- 
tion, we describe in Section |l| the details of the diffract- 
ing structures we consider, and the model we use to 
solve numerically the corresponding diffraction problem. 
Then we develop in section III the numerical resolution 
method, and discuss a series of examples. Finally, sec- 
tion |^ is devoted to the conclusions. More technical as- 
pects concerning the mathematical background and nu- 
merical details are discussed in three appendices at the 
end of this paper. 



II. MODELLING THE BIPERIODIC 
STRUCTURES 

We consider a system made of dielectric layers and 
biperiodic gratings of resistive conducting plates, ended 
by an infinitely conducting plane (or the vacuum), lo- 
cated at a height z — 0. The structure is globally in- 
variant under the discrete translations of period (a, b) 
which define the grating, namely translations of the form 
x — > x + ma, y — > y + nb, m,n G in the xOy plane. 
The structure is illuminated by an incident monochro- 
matic field of the form 



ft 



= ^ I e -i(vt- 



The geometry of the problem is displayed in Fig [l| We 
shall generically denote by and H w the electric and 
magnetic fields in the j-th layer Zj < z < Zj+i, with 
electric permittivity tf, we shall also use the superscript 
+ or — , according to whether the field propagates in the 
direction of positive or negative z. From now on, the 
configuration of Fig. [l] will be refered to as configuration 
I. 



*Thomson CSF-Optronique, rue Guynemer, BP 55, 78283 
Guyancourt cedex, France. 

^Birkbeck College, Physics Department, University of Lon- 
don, London WC1E, UK. Permanent address: CPT, CNRS- 
Luminy, Marseille. 



1 



7£ 




Dielectric 



Metallic 



Resistive plate 



FIG. 1. Global geometry of the structure: a grating of re- 
sistive plates between stacks of dielectric media, upon a per- 
fectly conducting plane. The incident field propagates in the 
direction of negative z. 

Alternatively, we shall also consider the same struc- 
ture, but we remove the infinitely conducting plane at 
z = 0. The latter configuration will be called configura- 
tion II. 



A. Floquet Modes 

Taking into account the global invariance of the prob- 
lem, it is natural to introduce the associated Floquet 
decompositions. Let k = [k x , k y , k z ) be the incident 
wavevector. In a medium of permittivity e, let us set for 
all integers m, n 



oto = k x ^fl 
Po = ky^/e . 
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where a and b are the grating periods. For m, n integers, 
we introduce the corresponding Floquet modes 



£mn( X > V> Z ) ~~ 
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and the "planar" modes 
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where ~&mn an d ~^ mn denote the (complex vector) co- 
efficients of the expansion. The sum over m, n runs the- 
oretically from — oo to oo. In practice it has to be trun- 
cated to a finite index [—M,M] x {—N,N}. We now re- 
strict ourselves to the tangential electric and magnetic 
fields. We directly obtain from Maxwell's equations that 
the following matrix relations hold 



* mn mn 
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Here we have introduced the following 2x2 matrices: 
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In the following we set £ mn = 

Alternatively, we shall also make use of the expansions 
with respect to the planar modes <f> mn (x, y) in (^), which 
leads to the coupled waves, defined by 
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The propagation of such modes within the corresponding 
layer is diagonal, and wc have in particular 



~^mn( z 3 + l) = exp{±ij mn {z j+1 - Zj)}E!n* (Zj) . (10) 

Matching boundary conditions at a dielectric-dielectric 
interface is an easy task, since Floquet modes with dif- 
ferent indices are not coupled. Given one such interface 
between two dielectric media labeled by j,j + 1, at a 
height z = Zj, and equating the tangential components 
of the electric and magnetic fields, we obtain: 
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where for the sake of simplicity we have suppressed the 
explicit dependence on the height z = Zj. The matrix 
elements (which are themselves 2x2 matrices) c, c' are 
given by 



which form an orthonormal basis of the space of biperi- 
odic functions on the plane, with period (a, b). Then it is 
well known that such function satisfy Helmholtz's equa- 
tion, and that both the electric and the magnetic fields 
may be decomposed into those Floquet modes (see e.g. 
[p"l|). Therefore, we write in the j-th layer 

^%,y,z)^jt^^E+ n (x,y,z) +-t^S~ n (x,y,zfj, (4) 

m.n 
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Alternatively, we shall make use of the following R- 
matrices, which read 
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and the connection between the two formulations is given 

by | 
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B. Surface Elements and Boundary Conditions on 
the Conducting Plates 

Let us now describe the surface currents living on the 
conducting plates. We may expand such currents into 
Floquet modes 



y(x, y) = ^ ~Jmn4>mn(x, y) 
m , n 

and impose the boundary conditions. 



(16) 




FIG. 2. The grating of conducting plates, in the particular 
case of rectangular plates. 

Several approaches have been proposed for imposing 
boundary conditions. Among these, the integral formu- 
lations (e.g. Galerkin methods) are generally considered 
the most stable. To implement the Galerkin method, we 
need to introduce a family of functions defined on the 
plates. Let %jj pq {x,y) be such a family. If zp denotes the 
height of the interface supporting the conducting plates, 
we then write, at a height z = zp 



(17) 



The boundary conditions rely on three sets of equa- 
tions. First, the continuity of the tangential electric fields 
at all interfaces 



^ (J+1 (x,y,Zj) = ^\x,y,Zj) 



3 J > 
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allows one to connect the global electric fields on each 
side of the interface. Second, the discontinuity condition 
for the tangential magnetic fields: 



^ (j+1 \x,y,Zj) ~ li^ ) {x,y,Zj) = X7f(x,y,Zj) , (19) 



explicitely 

(20) 

Finally, the impedance boundary conditions, which read 
at a height z = Zp\ 

~3 {P+1) (x, y, zp) = ~3 {P) (x, y, z P ) = Z7?(x, y, z P ) , 

(21) 

(where j vanishes outside the conducting plates) require 
a special treatment. It has been observed by several au- 
thors that such conditions cannot be imposed pointwise, 
because this leads to unstable systems. Several alter- 
natives have been proposed and tested (see for example 
||). The most stable solutions rely on the use of in- 
tegral formulations, obtained by considering either line 
integrals of the above equation, or a Galerkin formula- 
tion. We limit ourselves to the latter, which leads to a 
finite number of integral equations, obtained by testing 
Eq. (|2l| ) against suitably chosen basis functions tp pq (x, y) 
(see Appendix ^ for some possible choices) . 



C. The Coupled System 

Let us start with the case of configuration I. Taking 
into account the above remarks, we are led to the fol- 

lowing formulation. We denote by E and E the in- 
cident and reflected electric fields respectively, and we 
recall that we have denoted by P the index of the in- 
terface containing the plates. In order to avoid as much 
as possible numerical problems, we limit ourselves to a 
formulation involving the so-called i?-matrix propagation 
formalism j||5| (see Appendix |b] for a short account of 
the method). 

Using the i?-matrix propagation scheme, we can obtain 
R matrices for the stacks below and above z ~ zp. For 
example, we obtain a relation of the form 
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e! r 

mn 
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where the 2x2 matrices T and R are the stack equivalent 
transmission and reflection matrices respectively. Simi- 
larly, the i?-matrix algorithm below the grating of plates 
yields a matrix relation of the form 



(23) 
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The remarkable point with such a formulation is that 
it only involves small matrices, since modes with differ- 
ent indices m, n are not coupled. The only place where 
coupling between Floquet modes occurs is at a height 

z = Zp. 

The case of configuration II requires only minor mod- 
ifications. Eq. (^2|) is still valid. For the stack below the 
grating of conducting plates, we have to replace Eq. (|23| ) 
with 



(P)+ 



T'++ R 
R'+ T' 




(25) 



where mn are the Floquet coefficients of the transmit- 
ted field. Therefore, Eq. p^ ) is to be replaced with 

^mn = =- /V "~ 1 i?Lrl ■ (26) 

The rest of the formalism is unchanged. 



III. RESOLUTION AND NUMERICAL RESULTS 

A. Resolution of the Coupled System 

We now consider the practical resolution of the sys- 
tem we have obtained above. We consider approxima- 
tions of the fields with (2N + 1)(2M + 1) Floquet modes 
i? mn , m = —M, . . . M, n = —N, ... TV, and approxima- 
tions of the currents with PQ surface elements ip . The 
boundary conditions lead to three systems of equations 

involving the three sets of unknowns: ^^ + , ~E* m n ^ + 

—± ... -pt( p )+ 

and j mn - Eliminating E mn , we first obtain 

^mn ~=(1+N) 1 ((l + R-mn) ^Ln '+ T m«^m 



Inserting this result into (|f9|), we get 



(27) 
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where we have set 
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Eventually, we are led to a system of the form 

U mn ^^ mn \ r riinO mn j (31) 

where the 2x2 matrices U mn and V m „ are defined by: 



U mn — T, mji + (l + R^ m ) A mn B m „ , (32) 
V mn = Z-(l + R+;)A- 1 n X. (33) 

The system © is to be solved numerically, using a 
Galerkin procedure. Let ip pq (x,y) be a basis of func- 
tions defined on the plate, with appropriate boundary 
conditions. Using the expansion (|17|), we get 



~3 mn — ^ ] jpq 1p V q\mn i 



(34) 



where 

V* pq-mn = pqi ^mn) 



i> p q (x,y)(p mn (x,y)dxdy , 
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and where the star "*" denotes complex conjugation. 
Taking the scalar products of equations (|3f|) with the ba- 
sis functions ip pq (x, y), we obtain a system of the form 



^Pq ^ ] Vpq;p'q' jp'q' 
p',q' 



(36) 



where U is a vector of length PQ and V is a PQ x PQ 
matrix given by 



U vq — S ' (XJ mn E n 



V' pq-mn ) (37) 
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Eq. ( J36| ) is solved numerically (more details are given in 
Appendix ^|). Once the current J is known, one recovers 

directly the fields ~E~ m n ^ + using Eq. ( ^8| ) and then the 
reflected field ~E , from Eq. d2S 



B. Numerical Results 

Our main goal is to exhibit absorption bands, and to 
analyze the influence of some specific parameters on the 
location of the maximal absorption. More precisely, we 
focus on the influence of the resistive impedance Z and 
the ratio size of resistive plates/period. In addition, we 
show that the location of the absorption band essentially 
does not depend on the incidence angle. We work with a 
TM polarization for the incident field (in fact the results 
are weakly dependent on the polarization) . 

We consider a series of configurations, in which we vary 
individually these parameters, in the frequency domain 
IGHz — lOGHz. In all the figures, we plot the reflectivity 
(i.e. the ratio of reflected flux by incident flux) as a 
function of the incident frequency, and in the case of 
configuration II we also plot the transmittivity (i.e. the 
ratio of transmitted flux by incident flux) . 
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We start with the case of square resistive plates of side- 
length a', with variable impedance. The period of the 
grating is set to a = 10mm in both the x and y directions. 
The grating is supported by a dielectric stack of height 
z = Amm and complex refractive index e = 10 + 2i, it- 
self supported by an infinitely conducting plane (a simple 
case of configuration 7). Since the resistive plates are in 
that case square plates, we use Fourier-type decomposi- 
tions as described in Appendix A 1 for the decomposition 
of the surface current. The numerical results displayed 
below have been obtained using 17 x 17 Floquet modes 
and the same number of Galerkin modes. 
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Frequency (GHz) 

FIG. 3. Square plates in configuration I. Reflectivity as a 
function of the incident frequency, for various values of the 
plate's size: a' — 1,2, 4, 7, 9mm. The period a is kept fixed 
to a — 10mm. 6 — and Z = 0. 

We show in Figure || the reflectivity as a function of 
the incident frequency, for several values of the ratio a' /a. 
The computed values are indicated with symbols, and in- 
termediate values have been obtained using cubic spline 
interpolation. In all cases, a significant absorption band 
is observed. In addition, the critical frequency (i.e. the 
frequency at which reflectivity attains its minimum) de- 
creases as the ratio a' /a increases, and the width of the 
absorption band narrows. 

In the considered case, the plates are perfectly con- 
ducting. We nevertheless observe a strong absorption in 
a specific frequency range. Such a phenomenon is gener- 
ally coupled with the excitation of a leaky surface wave. 
The surface wave may be given an interpretation in terms 
of complex poles or zeroes of a scattering matrix (see B 
for details on the scattering matrix, and || for an anal- 
ysis of the role of zeroes and poles). The poles of the 
scattering matrix give the propagation constant of the 
leaky waves, which propagate along the surface of the 
biperiodic grating. The leaky wave is evanescent, as its 
energy decreases in the direction normal to the surface of 
the structure. The imaginary part of the pole gives the 
damping of the wave. The excitation of a leaky wave is 
a resonance phenomenon at a particular frequency. Fig- 
ure [| shows a spectacular phenomenon. A highly reflect- 



ing capacitive grating (with an important ratio a' /a) can 
absorb an incident plane wave in totality. Thanks to this 
absorption by a leaky surface wave propagating along the 
grating, we can control the absorption band of a classical 
Dahlenbach absorber layer which consists of a thick ho- 
mogeneous lossy layer backed by a metallic plate. When 
the ratio a' / a tends to zero we have checked that the min- 
imum of reflectivity is obtained for the same frequency 
as in Figure ||. To adjust the absorption band of the ab- 
sorber, we can deposit a biperiodic capacitive reflecting 
grating on the Dahlenbach layer. By doing so, we com- 
bine the properties of the biperiodic grating with those 
of the lossy layer. Notice in particular that it is possi- 
ble to decrease the thickness of the layer by adding such 
a biperiodic structure, to obtain the critical absorption 
frequency of the initial Dahlenbach structure. 



Square plates 




Frequency (GHz) 

FIG. 4. Square plates in configuration I. Reflectivity as a 
function of the incident frequency, for various values of the 
plate's impedance: Z = 10, 30, 100S1/D. The period a and 
the plate's size a' are kept fixed a — 10mm and a' = 7mm. 
6 = 0. 

We show in Figure [| the reflectivity as a function of 
the frequency of the incident beam, for several values of 
the impedance Z . The configuration corresponds to the 
case of Fig. ^ with a' = 7mm, and a significant minimum 
in the reflectivity is observed for a certain value of the 
frequency. This critical value is seen to be an increasing 
function of the impedance of the conducting plates. 

In Figure ||, the patches of the grating are not perfectly 
conducting any more. In that case, the absorption fre- 
quency and the bandwidth increase with the resistivity of 
the patches. To obtain a required absorption band, it is 
therefore possible to combine the effects of the geometry 
(here the ratio a 1 /a) and the effect of the conductivity. 
This provides extra flexibility to the filter design. 
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Square plates 




Frequency (GHz) 

FIG. 5. Square plates in configuration I. Reflectivity as a 
function of the incident frequency, for various values of inci- 
dence angle: 8 = 10, 30, 60 deg. The period a and the plate's 
size a' are kept fixed a = 10mm and a' = 7mm. Z = 0. 

We show in Figure || the reflectivity as a function of the 
frequency of the incident beam, for several values of the 
incidence angle 6, for the same configuration as before, 
i.e. a configuration exhibiting a well defined absorption 
band. These results (and other tests of intermediate in- 
cidence angles, are not reproduced here to simplify the 
plot) show that the critical frequency value depends very 
weakly on the incidence angle (at least for angles smaller 
than 45 deg). 

The same computations have been performed with re- 
sistive plates of various shapes. We display here the re- 
sults obtained when the square resistive plates in Fig. [| 
are replaced with cross-shaped ones, of the same size. By 
this we mean that the crosses lie within a square of the 
sidelength a', and are made of five identical squares of 
sidelength a'/3. For this cas e, w e used the surface ele- 
ments described in Appendix A 2 , and as before we take 



17 x 17 Floquet modes, and the same number of Galerkin 
modes. 

The numerical results, displayed in Figures § and 
show a similar behavior to the previous case: a well de- 
fined absorption band is clearly seen, and the critical 
frequency again depends on the ratio a' /a and on the 
impedance Z. Again, the location of the absorption band 
depends only weakly on the incidence angle (the numer- 
ical results, not given here, are very similar to those dis- 
played in Fig. ||). The only significant difference which 
may be observed is a broadening of the absorption band 
in the case of cross-shaped plates, and a second minimum 
occurs for large a'. 
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FIG. 6. Cross-shaped plates in configuration I. Reflectivity 
as a function of the incident frequency, for various values of 
the plate's size: a' = 1,2,4,7,9mm. The period a is kept 
fixed to a = 10mm. 9 = and Z = 0. 
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FIG. 7. Cross-shaped plates in configuration I. Reflectivity 
as a function of the incident frequency, for various values of 
the plate's impedance: Z = 10, 30, lOOfi/D. The period a and 
the plate's size a' are kept fixed a = 10mm and a' = 7mm. 
6 = 0. 



Similar computations have been made with configura- 
tion II. We display in Fig. || (reflexion) and Fig. ^| (trans- 
mission) the results obtained with systems identical to 
those considered in Figures ||-|5|. We observe that in such 
a configuration, the reflexion is small and almost constant 
above bGhz, it increases slightly with a' . The transmit- 
tivity shows maximums at frequencies corresponding to 
the minimums in configuration I. 

These two figures show the importance of the conduct- 
ing plane at z — 0. The well-defined absorption band 
appears only in that case. The excitation of the leaky 
wave and the corresponding absorption occurs only for 
structures ended by a conducting plane. 
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FIG. 8. Reflexion for square plates in configuration II. Re- 
flectivity as a function of the incident frequency, for various 
values of the plate's size: a' = 1,2,4,7,9mm. The period a 
is kept fixed to a = 10mm. = and Z = 0. 




Frequency (GHz) 

FIG. 10. Reflexion for square plates in configuration I. Re- 
flectivity as a function of the incident frequency, for various 
values of the plate's size: a' = 1,2,4,7,9mm. The period a 
is kept fixed to a — 10mm. 9 = and Z = 0. 
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FIG. 9. Transmission for squire plates in configuration II. FIG - 1L ^flexion for square plates in configura- 

Reflectivity as a function of the incident frequency, for various hon 1 Reflectivity as a function of the incident fre- 
values of the plate's size: a 1 = 1,2,4,7,9mm. The period a ^ nenc ^ for various values of the resistive impedance Z: 
is kept fixed to a = 10mm. 6 = and Z = 0. Z = 10, 25, 50, 75, lOOfi/D. The period a is kept fixed to 

a = 10mm. a = 7mm and 6 = 0. 



Next, we consider a second system (in configuration I), 
in which the resistive plates are located upon a double 
layer of dielectrics. The first dielectric (upon which the 
plates are located) has electric permittivity e = 5, and 
the second layer has electric permittivity e = lb -\- il&o / v 
with a frequency dependent imaginary part. Here the 
constant a is set to a = 10s /m, and the frequency v is 
expressed in GHz. 

The results are displayed in Figures [10] and [Uj As be- 
fore, an absorption band is clearly seen on Figure [l^, 
when a! is above 4mm, whose critical frequency and 
bandwidth decrease as the sidelength of the plates in- 
creases. In addition, for small plates, the reflectivity has 
a constant behavior close to zero above 8GHz. Figure |TT| 
shows that in such a configuration, the critical frequency 
depends weakly on the value of the impedance, but the 
bandwidth is an increasing function of the impedance. 



IV. CONCLUSIONS AND PERSPECTIVES 

We have studied and described a series of configura- 
tions involving dielectric stacks and arrays or resistive 
plates which produce well-defined absorption bands, with 
controllable absorption frequency. The critical frequency 
has been shown to be strongly influenced by the ratio 
period/plate-size, which therefore provides a good con- 
trol parameter. The impedance of the resistive plates has 
been shown to allow the control of the critical frequency. 

Our approach is based on a Floquet (or Rayleigh) de- 
velopment of the electromagnetic fields within the differ- 
ent layers of the structure, and a Galerkin approximation 
of the surface currents. Multilayers more complex than 
the ones we considered here may be described by the 
formalism of this paper as well. 
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In light of the numerical experiments we have per- 
formed, it is possible to combine the different parameters 
(namely the ratio a' /a, the geometry of the patches and 
the conductivity of the patch material) to obtain opti- 
mized absorbing structures from a quite standard biperi- 
odic grating. The use of absorption by a leaky surface 
wave can improve a classical Dahlenbach structure. 
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APPENDIX A: THE SURFACE ELEMENTS 

Depending on the geometry of the conducting plates, 
several different bases of surface elements may be used. 
In all cases, the finite number of basis functions we are 
forced to consider limits the precision of the approxima- 
tion of the current. 
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1. Rectangular Plates 

To start with, we consider the case of rectangular 
plates, as shown in Fig. || above. In such cases, the best 
choice for surface elements is provided by a Fourier basis: 
we set 

V> P9 V) = sm ^rF + l a \ cos — [y + ±b }e x 

+^cos— [# + Aa']sin— [y+|6']ej, , (Al) 

^ q M (x,y) = y sin-^- [x+\a'}cos— [y + \b']e x 

_pw co m^ x + i,j gin + i fo ']g* , ( A2 ) 

a' a' z b' 
Therefore, the Floquet modes of the surface current may 
be written as 

~3 1 mn — (jpq ^ pq.mn + 3pq ^ pq,mn , (A3) 

and the scalar products 

~$pq,mn = (^f , <f>mn) , (A4) 
V* pq,mn = {^pq i 4>ran) i (A5) 

may be computed analytically. 

For other special geometries, such as disks or elongated 
disks, it is possible to design appropriate basis functions 
to describe the current density on the resistive plates (in 
the case of disks, such basis functions are linear combi- 
nations of Bessel functions). However, it is also desirable 
to have basis functions which can describe arbitrary ge- 
ometries. This is the purpose of the surface elements 
described in the next subsection. 



8 



2. Arbitrary Plates 

For conducting plates with arbitrary geometry, we are 
forced to use "all purpose" basis functions, which we shall 
call surface elements. Such basis functions have been 
considered by several authors under the name of rooftop 
functions. It follows from the analysis in || that rooftop 
functions often provide faster and better conditioned nu- 
merical schemes than classical alternatives (the so-called 
surface-patch and triangular patch functions). The first 
step for the construction of such surface elements is a 
discretization of the plate. For the sake of simplicity, we 
restrict to a uniform square discretization, with period r. 
Consider the characteristic function 



1 if < X < T 

elsewhere 



and the Schauder function 

1 + f if - t < z < 
A(as) = < 1 - 1 if < x < t 







elsewhere 



Then set 



V'pqOz, V) = X(x - pr)A(y ~ qr) , 

^p q (x, y) = A(z - pr)x{y - qr) , 
and finally 



^ pq (x,y) 



pq e y 



(A6) 

(A7) 

(A8) 
(A9) 

(A10) 



The surface elements we consider will be those functions 
tpp q (x,y) and ip y q (x,y) such that their support is com- 
pletely included in the support of the plate. Clearly, the 
smaller r the better is the approximation of the current, 
but the higher the complexity of the numerical problem. 



APPENDIX B: R-MATRIX PROPAGATION 

We describe briefly the i?-matrix propagation scheme 
as we used it in our simulations. Clearly, the simplest 
approach amounts to consider the direct product of the 
C matrices given in Eq. (|ll]), which yields directly a C 
matrix for the whole structure. As stressed by various 
authors, such a scheme turns out to become rapidly un- 
stable as the depth of the structure grows. 



1 



= R 



1 



T 



E%) E© + ( Zj ) 



Let us consider a multilayered medium with interfaces 
at heights z p , . . . , zf, and assume that we are given an 
interface i?-matrix of the form given in Eq. ([lj]). Then, 
one easily verifies that 



(i)+, 




where we have set 



zA 



— t Lj 5 t 


t Lj ; 


= r+-L* ; f"+ 


= r-+ , 


exp{ij mn (zj + i 











ae g+1)+ 







(Bl) 



(B2) 



FIG. 13. Illustration of the Ji-matrix propagation algo- 
rithm: the stack ii-matrix. 

Suppose now that we are given a stack i?-matrix for 
the stack [j + 1,F]: 




(B3) 



(j)± 



Zj-\) for the 



where we set by default E 
sake of simplicity. FromEqs. (Bl) and (B3), little algebra 
gives the expression of the coefficients of the stack matrix 
for the stack [j + 1,F]: 



^ R (z F ) = T++{l- R+-f 
+ (R~++T 



-+(1- R+-f-+)- 1 f-+T— )^( 



zf) 



Zi 



.i) = (f +- +t~(l- R+-f-+)- 1 f- 



£? 0) ( 



■3-D 



'(I- R 



+ -S-+' 1 -lr r 



zf) 



The above equations provide a simple iterative algorithm 
for computing the global R- matrix for the stacks [zp, Zf] 
and [zo,zp]. This algorithm is known as the R-matrix 
propagation algorithm, and has been analyzed by various 
authors. We refer to l0,S,[7lJTo|l for more details. 



APPENDIX C: NUMERICAL ASPECTS 



FIG. 12. Illustration of the J?-matrix propagation algo- 
rithm: the role of an interface .R-matrix. 



We give here more details on the numerical methods 
used to solve the complete problem. As stressed before, 
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most of the matrices used in the scheme are 2x2 ma- 
trices, which are easy to handle. In addition, the use of 
i?-matrix propagation algorithm prevents us from devel- 
oping numerical instabilities when computing products 
of such matrices. 

The main part of CPU is used for solving Eq. (^6|). Sev- 
eral methods have been tested for that problem (which 
has also been studied by various authors). The numer- 
ical results presented here have been obtained by using 
an inversion method based on L/7-decomposition, with 
left and rigth equilibrations of the matrix. A fortran im- 
plementation of such a method is available in the LA- 
PACK library (see ||). Alternative methods may be 
found in the literature, such as (complex) biconjugate 
gradient methods or FFT-based methods. 
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